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Abstract 

The dependence of ionospheric 0 + escape flux on electromagnetic energy flux and 
electron precipitation into the ionosphere is derived for a hypothetical ambipolar 
pick-up process, powered the relative motion of plasmas and neutral upper atmos- 
phere, and by electron precipitation, at heights where the ions are magnetized but 
influenced by photo-ionization, collisions with gas atoms, ambipolar and centrifugal 
acceleration. Ion pick-up by the convection electric field produces "ring-beam" or 
toroidal velocity distributions, as inferred from direct plasma measurements, from 
observations of the associated waves, and from the spectra of incoherent radar ech- 
oes. Ring-beams are unstable to plasma wave growth, resulting in rapid relaxation 
via transverse velocity diffusion, into transversely accelerated ion populations. Ion 
escape is substantially facilitated by the ambipolar potential, but is only weakly af- 
fected by centrifugal acceleration. If, as cited simulations suggest, ion ring beams re- 
lax into non-thermal velocity distributions with characteristic speed equal to the lo- 
cal ion-neutral flow speed, a generalized "Jeans escape" calculation shows that the 
escape flux of ionospheric 0 + increases with Poynting flux and with precipitating 
electron density in rough agreement with observations. 
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1 Problem 

2 Heating and ablation of ionospheric plasma by solar wind energy is an important 

3 atmospheric loss process that shapes Earth’s magnetosphere during space storms, 

4 adding substantial plasma pressure to the magnetosphere [Moore and Horwitz, 

5 2007]. The rate of removal of the atmosphere is non-threatening over human time 

6 scales, but is representative of a widely applicable space plasma process that played 

7 a role in removing much of the atmosphere of Mars. One possible agent of such heat- 

8 ing and acceleration is plasma waves driven by magnetic field-aligned electric cur- 

9 rents that transmit stresses to the auroral ionosphere. The heating should then be 

10 maximal in the auroral current sheets at about 1 Re altitude. High altitude heating is 

11 effective in raising the flow speed of ion outflows. However, the escaping mass flux 

12 is determined by energy inputs at lower altitudes [Leer and Holzer, 1980; Moore 

13 and Horwitz, 2007] in the topside F region ionosphere, about 250-1000 km altitude. 

14 Strangeway et al. [2005] and Zheng et al. [2005] found that two local magneto- 

15 spheric factors are well correlated empirically with outflow flux: i] the precipitating 

16 magnetospheric electron density and ii] the DC electromagnetic [or Poynting] flux 

17 into the F region ionosphere. Plasma outflow resulting from soft electron precipita- 

18 tion can be understood [Caton et al. 1996] as the effect of ambipolar coupling of ions 
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1 and electrons. Heavy ion plasmas are also driven out of the ionosphere by wave en- 

2 ergy dissipation and pressure acquired from the heating [Zeng and Horwitz, 2007]. 

3 The source or sources of effective waves has remained elusive, but the Strangeway 

4 and Zheng results cited above strongly suggest that the waves may be driven by 

5 ionospheric convective motions through the neutrals, since that is the principal con- 

6 sequence of DC Poynting flux into the ionosphere. 

7 The requirements for ion escape flows are demanding and have been widely 

8 thought to preclude significant escape of heavy ions owing to Joule or frictional en- 

9 ergy dissipation. The gravitational escape velocity of 0 + being approximately 11 

10 km/s in the F region, bulk escape of ionospheric 0 + is thought to require the average 

11 0 + ion to be given over 10 electron volts, implying a temperature of ~100,000 K. Lo- 

12 ranc and St-Maurice [1993] studied transient frictional heating as a method of pro- 

13 ducing ion up-flows, and they showed that transient effects produce non- thermal 

14 velocity distributions that are not well described in fluid approaches. However, they 

15 did not concern themselves with processes that turn up-flows into outflows by ac- 

16 celerating ions to escape velocity. 

17 Ions are significantly accelerated transverse to the local magnetic field, suggesting 

18 that a wave generation and dissipation mechanism must operate to raise the gyra- 

19 tion velocities of ionospheric ions. Ionospheric convective motions driven by the ac- 
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1 tive magnetosphere are typically in the range of 1-3 km/s in active auroral situa- 

2 tions. This falls substantially short of the nominal gravitational escape velocity, 11.2 

3 km/s, but significantly raises ion scale heights. To get outflows, modern ionospheric 

4 simulations [e.g. Barakat et al. 2003, 2006] usually invoke a combination of electron 

5 precipitation or heat flow, with application of empirically determined energy inputs 

6 often characterized as "wave particle interactions", for example by Zeng and Hor- 

7 witz [2007]. The problem addressed here is the energy source and production proc- 

8 ess of the waves and acceleration required to produce observed outflow fluxes. 

9 Strangeway et al. [2005] and coworkers cited therein also identified Alfven waves, 

10 generated in the magnetosphere and dissipated in the ionosphere, as agents of 

1 1 plasma outflow, and showed that outflow flux is well correlated with the Poynting 

12 flux of AC or wave energy propagating into the ionosphere along magnetic flux 

13 tubes. Such waves may exist in the ion resonant frequency range, or may non- 

14 resonantly power outflows via the ponderomotive effect [Khazanov et al. 1998; Gug- 

15 lielmi and Lundin, 2001]. As important as these effects may well be, they remain 

16 problematic to derive from global magnetospheric simulations, owing to the inher- 

17 ent time resolution or other limitations of such models. In the future, we can antici- 

18 pate a capability to specify magnetospheric wave spectra based on such models that 

19 will expand from DC to higher frequencies, and will support making such a connec- 
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tion with ionospheric dissipation and outflow. For the present, however, we limit 
ourselves in this paper to DC effects that can be extracted from current global simu- 
lations of the magnetosphere-ionosphere system. 

Inasmuch as simple escape velocity arguments have been invoked to argue against 
important outflow effects of ionospheric convection and frictional heating, we adopt 
the same approach, to see if such arguments can be rebutted via a hypothesis based 
on ion pick-up by the convection electric field. That is, we use a generalized "Jeans’ 
escape" calculation for ions in which the gravitational potential is reduced by the 
ambipolar potential and by a centrifugal potential, while the non-thermal ion veloc- 
ity distribution of pick-up ions (PUI) is assessed and used to compute the fraction of 
the ions that are freed from gravity as reduced by this combination of effects. Sub- 
sequent sections deal with the ambipolar and centrifugal potentials, the convective 
pick-up process, convective transfer of energy, pick-up ion ring beam relaxation, 
generalized Jeans escape, and the calculation of 0 + escape flux. A discussion section 
is followed by conclusions. 

Ambipolar Potential 

The ambipolar electric field couples the energy of electrons, which would otherwise 
readily escape gravity, to the heavy and slow ions. The ambipolar potential simulta- 
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neously traps the electrons and reduces the gravitational binding of ions, allowing a 
fraction of them to escape. In steady state, charge balance adjusts the total ambipo- 
lar potential drop such that the net escape fluxes of electrons and ions to the magne- 
tosphere are equal. When a superthermal population of electrons has energies sub- 
stantially larger than the ion gravitational binding energy, their fractional density is 
effective in controlling ion escape, as is ion heating. Khazanov et al., [1997] consid- 
ered the dependence of the ambipolar field and its integrated potential on the pres- 
ence of un-thermalized photoelectrons. They found that the photoelectron source 
produces fractional densities of superthermal electrons in the range of 0.02 to 0.04 
%, corresponding to a total ambipolar potential above 500 km of ~ 4-5 V, and this is 
consistent with observed polar wind H+ escape with relatively little 0 + escape [Abe, 
et al., 1993]. The photoelectrons were assumed to have a characteristic energy of 20 
eV, but the results were only weakly influenced if the mean photoelectron energy 
was increased to 40 or 60 eV. As the fractional density of photoelectrons was varied 
from 0.01% to 1%, the ambipolar potential varied from about 2.9 V to 7.3 V, the 
variation appearing as a transition between asymptotes, as shown in Figure 1. The 
lower asymptote corresponds to the ambipolar potential for ion escape to match 


thermal electron escape. The upper asymptote corresponds to the ambipolar poten- 
tial required to free enough 0 + ions to match the escape of the superthermal elec- 
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trons. Because the superthermal electrons are assumed to have characteristic en- 
ergy significantly larger than the 0+ gravitational escape energy, the required ambi- 
polar potential is insensitive to superthermal electron characteristic energy. Be- 
cause ambient 0 + ions have thermal energy substantially less than their gravita- 
tional escape energy, the required ambipolar potential saturates as it approaches 
bulk thermal 0 + escape, which increases steeply with incremental further increase 
in ambipolar potential. Enhanced heavy ion escape owing to other factors such as 
ion heating or acceleration effects will similarly limit or suppress the ambipolar po- 
tential to a value less than that otherwise needed to balance ion and superthermal 
electron escape. 

Secondary electrons are produced by primary electrons of magnetospheric origin 
and have the same interaction with ionospheric outflows as photoelectrons, via the 
ambipolar potential. The shape of the low energy secondary electron distribution 
also steepens toward the low thermal energies because of collisional plasma ther- 
malization processes, but may differ in detail from the photoelectron spectrum. 
Nevertheless, because their characteristic energy is greater than the 0 + binding en- 
ergy, the ambipolar potential is insensitive to the spectrum details and is driven by 
the superthermal electron partial density, as shown in Figure 1. 
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For present purposes, we assume that the precipitating electron density identified 
by Strangeway et al. [2005] characterizes the auroral secondary electron population 
and is proportional to the fractional superthermal electron density of Khazanov et 
al., [1997], with the constant of proportionality regarded as a free parameter to be 
determined. Electron precipitation was evaluated for E > 50 eV in the former study, 
while the latter study considered the entire superthermal electron distribution, 
taken to have a characteristic energy of 20-60 eV. Since the density of electrons >50 
eV is necessarily much less than the full density of superthermal electrons without 
such a lower limit, the total density of electrons that produces the correct ratio 
Np/Ntot will be smaller than typical local electron densities by a factor that depends 
upon the actual shape of the superthermal electron velocity distribution. Figure 1 
displays a hyperbolic tangent fit to the relationship found by Khazanov et al., [1997], 
together with the fitting parameters and their fitting uncertainties. This ambipolar 
potential reduces the gravitational escape velocity for ionospheric ions, as a function 
of precipitating electron density, according to the relationship: 


where G is the gravitation constant, Me and Re are the mass and radius of Earth, H exo 
is the exobase height, q is the elementary charge, and Mi is the ion mass. The present 
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work would clearly benefit from more comprehensive precipitating electron obser- 
vations without a lower energy limit. This is an important target for future observa- 
tional research in this area, because a plentiful supply of superthermal secondary 
electrons will cause the ambipolar potential to rise to the level where thermal ion 
escape is incipient and requires only a small amount of additional energy to become 
a large bulk escape. 

Centrifugal Acceleration 

A prominent source of energy for ion acceleration is convection driven by magnetic 
coupling with magnetospheric hot plasma motions. Terrestrial F-region convection 
is much slower than solar wind or magnetospheric flow speeds, owing to the large 
ratio of magnetic field intensities. However, ionospheric convection is comparable 
to or exceeds neutral gas thermal speeds. It also increases steeply with altitude, 
such that higher altitude plasmas on the same flux tubes may approach the gravita- 
tional escape speed, independent of ambipolar effects. These parameters are illus- 
trated in Figure 2. Here the Poynting Flux at 4000 km altitude has been set to 1, 10, 
and 100 mW/m 2 so the system is driven across the effective range observed in the 
Strangeway [2005] study, for a typical value [8 S] of dayside auroral ionospheric 
[Strange way, private communication]. 
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The increase in the convection speed with altitude shown in Figure 2 gives rise to a 
centrifugal acceleration effect [Horwitz etal., 1994; Demars etal., 1996]. Plasma 
ions flowing along rapidly convecting magnetospheric flux tubes are accelerated as 
both perpendicular and parallel energy are acquired from the convection electric 
field. Ions gain equal amounts of perpendicular and parallel energy, equal to the 
convection speed at any altitude, adjusted for gravitational accelerations parallel to 
the magnetic field [Horwitz et al., 1994]. Centrifugal forcing depends weakly on con- 
vection path and field curvature details, but for convection across the pole of a di- 
pole field (their equation 4b], at altitudes lower than a few Re it is approximately: 


F = 1 5 mV 2 — 

1 1 rrl V co nv 3 

Kr 


( 2 ] 


This centrifugal force can be expressed as the gradient of a potential: 


<D 


1 2 f r ] 

2 


( 3 ] 


This reduces the escape speed everywhere via the subtraction of an additional term 
from equation [1], yielding equation [4]. 
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The escape speed falls to zero at the levels shown in Figure 2, for different values of 
the Poynting flux (or convection speed), with zero ambipolar potential. Non-zero 
values of ambipolar potential clearly reduce the escape speed further. Above the 
convection (and ambipolar) dependent altitude where escape speed drops to zero, 
which we term the "centrifugal horizon," even cold ions continue to move upward 
and are accelerated further. This height is lowered by strong convection. However, 
in the topside ionosphere, below the exobase, the escape speed is not reduced by a 
substantial amount, and additional acceleration is needed to overcome gravitational 
binding, even when convection is very strong. Thus, as shown by Demars et al., 
[1996], centrifugal acceleration alone is not a significant effect on outflow flux, 
thought it is important in accelerating the outflows, with obvious effects on their 
circulation in the magnetosphere. It is included here mainly for completeness. 

Convective Pick-up 

The generation of ion-resonant diffusive waves is familiar from other situations 
where photo-ionization creates new ions with initial velocities of the parent cold 
neutral atoms, which are then "picked up" by the local plasma motional electric field. 
Pick-up is distinct from dissipation of horizontal currents in the lower parts of the 
ionosphere where ions are collision dominated, usually known as Joule or frictional 
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1 heating. The pick-up process operates at altitudes where collisions with atoms or 

2 molecules are rare compared with gyrations, and photo-ionization or charge ex- 

3 change dominate ion-neutral interaction, that is, above ~125 km altitude [Kelley 

4 and Heelis, 1989], a boundary to which we refer as the "isotropause." Above that 

5 level, magnetic anisotropies can be created and persist without being isotropized by 

6 collisions. Most of the F-layer is above the isotropause, including its "topside" ex- 

7 tending to the exobase, defined as the altitude where the mean free path exceeds the 

8 scale height of the scattering neutral gas. The exobase lies in the range from 500 to 

9 1300 km altitude [Banks and Kockarts, 1973], increasing linearly with the neutral 

10 thermospheric temperature during periods of active Joule heating. 

11 Pick-up ion (PUI) distributions are produced by convection of ions through the neu- 

12 tral gas, when ions are created in the neutral gas frame by photo-ionization [if sun- 

13 lit] or by charge exchange between ions and neutrals, or to a lesser extent by scat- 

14 tering collisions of the drifting ions with stationary neutrals. Examples include the 

15 interstellar gas flow throughout the solar system [Bogdan et al., 1991], or the local- 

16 ized gas sources at un-magnetized planets, satellites, and comets [Szego et al., 2000]. 

17 PUI are the accelerated cold ions, after a gyro-period, during which they are picked 

18 up by the convection electric field as seen in the neutral frame, gaining drift and gy- 

19 ration speed each equal to the convection speed [Kivelson & Russell, 1995, p.206]. 
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The result is a ring or toroidal velocity distribution, or a "ring beam", rooted in the 
neutral frame and encircling the convecting plasma frame with a diameter of twice 
the plasma convection speed in velocity space. This mechanism was studied in an 
ionospheric context during the 1970‘s, reviewed by St-Maurice and Schunk, [1979], 
revisited by Barakat et al., [1983], Winkler et al. [1992], Wilson [1994], Barghouthi 
et al., [1994], and Hubert and Leblanc [1997]. PUI are an inevitable consequence of 
moving driver plasma that is magnetically linked to a partially ionized "anchor" gas, 
in the altitude range where collisions are important but infrequent relative to the 
local gyro frequency. 

The geomagnetic field couples magnetospheric convection into the ionosphere as 
shown in Figure 2, but the plasma motion is opposed by drag owing to ion collisions 
with the gas atoms, which also demagnetizes the ions and allows transverse current 
to flow in the low F region, principally below the isotropause. Maxwell stresses 
(JxB) drive ionospheric plasma flow and retard the conjugate magnetospheric flow, 
transmitted by field-aligned electric currents. The ionospheric flow is powered by 
the Poynting energy flux from the driver plasma into the convecting ionospheric 
plasma. Thus, either Poynting flux (Fp) or plasma flow can be regarded as the driver, 
but we adopt Fp as the driver to facilitate comparison with the results of Strangeway 
[2005]. Given fixed Fp, ionospheric convection speed depends on the ionospheric 


13 


Moore & Khazanov - Mechanisms of Ionospheric Mass Escape 


1 conductance. But in either case, the ionospheric height integrated conductance (Ep) 

2 is the load on the electric circuit or on the driving motions, and the convection speed 

3 corresponding to a particular value of Fp flowing into the ionosphere is given by: 

4 v„,-j4i/VS 7 C5) 

5 Conductance distributions given by Ridley et al. [2004] are representative, indicat- 

6 ing a range of £p from 2-20 Siemens, representative of the entire dayside region and 

7 polar cap under a wide range of conditions, and compatible with 8 S as a typical 

8 value. For given Poynting flux, convection speeds will be correspondingly lower in 

9 the low latitude dayside or night time auroral zone, but Poynting flux may also be 

1 0 larger, with compensating effect on the magnitude of ionospheric convection speed. 

11 Ring Beam Relaxation 

12 In practice, ring beams have proven difficult to observe, consistent with rapid re- 

13 laxation into non-thermal but transversely accelerated distributions. Ion species 

14 self-collisions are expected to erode and thermalize ring beams [Barghouthi et al., 

15 1994] on collisional time scales. On faster plasma time scales, unstable velocity dis- 

16 tributions will be eroded away by plasma waves. Any velocity distribution with 

1 7 positive df(v) /dv has free energy available to drive wave growth that removes the 
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1 free energy feature by diffusing particles from higher to lower phase space densities 

2 and predominantly, though not exclusively, from higher to lower energy. The diffu- 

3 sion occurs along characteristics specific to the driven wave modes. The magnitude 

4 of df(v ) /dv and corresponding wave growth driven by the ionospheric PUI feature 

5 are compounded in the topside ionosphere by the highest phase space densities and 

6 lowest characteristic speeds anywhere in the space environment. As a result, one 

7 expects to observe distributions that are marginally stable with df(y)/dv < 0. St.- 

8 Maurice and Schunk [1979] cited plasma observations of the ring feature without a 

9 pronounced minimum at zero velocity. The ratio of the ring speed to neutral thermal 

10 speed (often referred to as D*) never exceeded 1.5 in that study, even when ion 

11 temperature suggested a D* of 2 to 3 or more. Space and time variability of the iono- 

12 sphere is also a significant problem for observations, making it difficult to capture a 

13 ring beam feature [Moore et al., 1996]. High altitude long-lived observations indicat- 

14 ing a core ring beam features have been reported [Moore et al., 1986], but are rare. 

15 The above facts are consistent with ring beam production of intense waves that 

16 quickly erase the velocity space instability responsible for them. 

17 In contrast, many observations have been made of "transversely accelerated ions" or 

18 TAI, and their close relatives, ion conics [Moore and Horwitz, 2007]. Conics develop 

19 as transversely accelerated ions move upward while slower ions remain relatively 
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stationary. The prevalence of auroral TAI and ion conics requires that transverse 
energy diffusion be faster than pitch angle diffusion in the auroral context. Figure 3 
shows observations by the Swedish Freja spacecraft [Andre etal., 1994]. The per- 
pendicular velocity distribution is a power law that could be approximated as expo- 
nential for low speeds. The thermal core of this distribution, if any, is at velocities 
lower than the range observed. At higher speeds, the distribution is well fit by a 
generalized Lorentzian or power law with (v/v s d) 2 - 2 as argument, as indicated in 
Figure 3 and identified in Table 1. Andre et al. [1994] simulated the effect of the ob- 
served ambient wave populations at Freja and concluded that the observed waves 
would heat ions so as to produce the observed hot tail distributions. The open ques- 
tion is that of the source of the waves that heat the ions, principally in the transverse 
direction. Local wave growth owing to a ring beam feature was derived by Post and 
Rosenbluth [1966], but in the context of the loss cone instability. Such waves are 
commonly observed in the topside ionosphere under conditions of strong convec- 
tion [St-Maurice and Schunk, 1979]. Significantly, a number of reports of incoherent 
scatter radar echo spectra have been interpreted as the result of backscatter from 
plasmas with a toroidal or ring beam velocity distribution [Suvanto et al., 1989; Kin- 
zelin and Hubert, 1988]. 
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Bogdan et al., [1991] interpreted shell distributions in interplanetary space as prod- 
ucts of the PUI process. Near comets, plasma analyzers encounter PUI populations 
[Neugebauer et al., 1990], and deeper in the coma, "thermalized" exponential or 
power law tails attributable to them [Mukai, 1986; Richardson et al. 1987]. Theo- 
retical work has treated the instability in the high velocity outer coma [Lee and 
Gary, 1991], and as a function of distance extending deep into the inner, collisional 
coma [Puhl et al,. 1993], using a velocity diffusion equation with empirically deter- 
mined coefficient in the latter case. However, the general self-consistent quasi-linear 
evolution problem has not evidently been solved. 

Quasi-linear plasma theory describes diffusion along characteristics that are circles 
centered on ± Va (Alfven speed] in the weak-field, high-speed, collisionless regime. 
Velocity diffusion along such circles centered near the origin in velocity space, rela- 
tive to the ion ring beam, takes the form of pitch angle diffusion and forms bi- 
spherical shell distributions [Szego et al., 2000], with little energy diffusion. Deeper 
in a cometary coma, ion flow velocities fall to much smaller than Va. Then, as for the 
Earth’s ionosphere, diffusion characteristics are circles centered at large parallel ve- 
locities, which cut transversely through the ring beam. Also, dispersive ion cyclotron 
waves may become important, with non-circular diffusion characteristics [Isenberg 
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1 and Vasquez, 2007]. In either case, velocity diffusion of the PUI feature is transverse 

2 to the magnetic field, and pitch angle diffusion is relatively weak. 

3 Puhl et al. [1993] simulated coma PUI ring beam "thermalization" and found a 

4 power law form that could also be approximated as exponential at low speeds. They 

5 found the "temperature" of the relaxed distribution to correspond to a thermal 

6 speed equal to double the local convection speed (their figure 6b], or the bulk plus 

7 gyration velocity of ions in the frame of the neutrals. However, this "temperature" 

8 was derived as a moment of highly non-thermal power law velocity distributions. 

9 We take the scale velocity for all distribution forms to be equal to the local convec- 

10 tion speed, representing the gyration energy in the convecting ion frame of refer- 

11 ence. The bulk convection speed is also used to evaluate the centrifugal potential. 

12 The above simulation results were supported by a relevant active experiment re- 

13 ported by Paterson and Frank [1989] and Gurnett et al. [1988]. They observed a 

14 power law tail of hot ions in the wake of the space shuttle, from the Plasma Diagnos- 

15 tic Package, and inferred that charge exchange between shuttle out-gassing and 

16 ionospheric ions produced new ions in the shuttle frame, which were then "picked 

17 up" into the ionospheric plasma frame, producing ring beams that generated the ob- 

18 served waves and subsequently turned into hot transverse ion distributions with 

19 measurable fluxes up to 100 eV. We find that the observed distributions are well fit 
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1 by a kappa distribution with kappa ~ 2.75 and a characteristic speed of ~ 3 km/s. 

2 These ion distributions remained transversely peaked and did not form spherical 

3 shells as is the case for cometary ion pick-up, consistent with the high Alfven speed 

4 in the ionosphere and diffusion principally transverse to the local magnetic field. 

5 From the above, we conclude that ion ring beam distributions relax through per- 

6 pendicular velocity diffusion, in the low speed, high-VA context appropriate to the 

7 ionosphere or inner coma of comets, both in theory and in practice. An extended 

8 transverse velocity distribution is expected to form, having power law dependence 

9 at higher speeds, with scale speed equal to the local convection speed. Incidentally, 

10 we note the ubiquitous presence of power law tails in space plasmas, which may be 

1 1 understandable if ion pick-up is similarly ubiquitous. In any case, we adopt this 

12 simulated and observed distribution to assess ion escape. We note that any other 

13 source of waves that produces power law tails scaled with the local convection 

14 speed would be equally effective in producing the results given below. 

15 Important assumptions of this approach should be noted: First, unstable waves are 

16 assumed effective in relaxing the distribution as soon as the ion convection speed 

17 becomes an appreciable fraction (10-30%) of the neutral thermal speed. This is the 

18 same condition required for a perceptible local minimum to form in the velocity dis- 

19 tribution owing to the gyration of newborn ions and thus seems justifiable. How- 
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ever, a full ring beam instability analysis must be done to substantiate this assump- 
tion. Second, it is assumed that the generated wave intensities, like other imposed 
wave intensities discussed above, have an insignificant ponderomotive effect on the 
ions [Khazanov et al., 1998]. Such a ponderomotive effect would enhance outflow 
above the values obtained here, so this is a conservative assumption that neverthe- 
less may cause the present results to underestimate outflow fluxes. Finally, we as- 
sume that the time scale for such power law distributions to relax to Maxwellian 
would be a collisional time scale much longer than the wave particle diffusion time 
scales responsible for creating the power laws. 

Generalized Jeans' Escape 

The standard assessment of 0 + escape is that such escape requires on the order of 
10 eV to be imparted to each escaping ion, so that it can overcome Earth’s gravita- 
tional binding. Equivalently, this requires each 0 + ion to be accelerated to 11.18 
km/s, the "escape velocity" for Earth. The typical temperature of oxygen in the iono- 
sphere being of order 3000 K, or about 0.03 eV, corresponding to a thermal speed of 
only about 300 m/s, typical oxygen ions must gain most of 11 km/s to escape. This 
argument is a direct application of Sir James Jeans [1904] approximate method for 
calculating the rate of loss of gas atoms from a planet. At the exobase, above which 
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the mean free path exceeds the scale height so that upward motion is free of colli- 
sions, the species density and thermal velocity, (that is, the species thermal or limit- 
ing flux) are multiplied by a factor corresponding to the fraction of the atoms with 
velocities exceeding the escape velocity. This simple method has been shown to ex- 
aggerate the H and He escape rates by about 3-30% [Brinkman, 1970], but this level 
of accuracy is more than sufficient for present purposes. 

We generalize the Jeans method, adapting it for ion escape in two ways. First, we 
consider non-thermal forms of velocity distribution in addition to thermal Maxwel- 
lians [see also Shizgal and Arkos, 1996]. Table 1 provides a listing of the forms used. 
Second, but just as important, we consider additional potentials that act on ions, 
namely the ambipolar electric field and centrifugal force, as discussed above. We 
estimate the escape flux as the product of the limiting flux for 0 + at the exobase and 
a factor corresponding to the fraction of the ions exceeding the exobase escape ve- 
locity. We treat velocity in one dimension, because perpendicular velocities become 
parallel velocities (in the source cone) for free ion motion above the exobase, con- 
serving the first invariant. 
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1 Oxygen Escape Flux 

2 We next use the "generalized Jeans escape" estimate described above to evaluate the 

3 dependence of ion escape on Poynting flux and electron precipitation, evaluating 

4 escape at the exobase. Also relevant is the 0-H density crossover, above which fast 

5 0 + charge exchange produces mainly cold H+ ions, rather than cold 0 + ions. This oc- 

6 curs near the exobase, with a similar dependence on thermospheric temperature 

7 [Moore, 1980]. Scattering collisions with neutrals continue to be important above 

8 this altitude, with the effect of maintaining and enhancing the PUI feature, as con- 

9 vection speed increases. Scattering collisions with neutrals interrupt ion drift mo- 

10 tions and randomize them about the neutral gas frame, such that they are again 

1 1 picked up by the electric field, albeit with less contrast between their gyrating drift 

12 motion and their random motions than for charge exchange. Scattering collisions 

13 with other ions, on the other hand, tend to relax the ring beams into transversely 

14 heated velocity distributions. We neglect these competing effects and indeed all 

15 such processes above an exobase taken at 500 km, regardless of Poynting flux. This 

16 conservative approach ignores the rise in the exobase height with thermospheric 

17 temperature, even as most of the Poynting flux goes into frictional heating of the 

18 neutral gas. It will therefore tend to produce underestimates of the escape flux dur- 
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ing very disturbed conditions that raise the effective exobase as the thermosphere 
heats up. 

The escape flux of 0 + is limited by photo-ionization sources and friction with ther- 
mospheric gas. The limiting flux has been computed via a full ionospheric topside 
simulation [e.g., Barakat et al., 1987]. It depends on many factors including solar ac- 
tivity, F-peak density, influenced both by photo-ionization and horizontal transport, 
and the topside scale height of the neutral gas, influenced by Joule heating. Both 
photo-ionization and charge exchange contribute to the production of ions from cold 
neutrals. Active solar conditions produce nearly 10 3 0 + cnr 3 s 1 at 125 km altitude 
[Banks and Kockarts, 1973, p.B-114]. Integration over a scale height of 100 km 
above that supports a flux of 1 x 10 10 cnr 2 s - 1 . A similar cold ion creation rate applies 
to 0 + -0 charge exchange [Banks and Kockarts, 1973, p.A-224], if we take no+ ~ 10 s 
cm 3 , and no ~ 10 7 cnr 3 [Kelley and Heelis, 1989, p.6] with 100 km scale height. We 
thus have an adequate source of cold ions for pick-up to provide an 0 + limiting flux 
(at 4000 km altitude) of ~ 2 x 10 10 cnr 2 s _1 , which is the largest flux observed by the 
Polar TIDE investigation [unpublished white paper]. This is also the lowest value 
consistent with the lack of apparent saturation evident in the Strangeway et al. 
[2005] results, shown as the background of Figures 4-6. Higher production would 
follow from the enhanced supra-auroral 0 densities derived from recent spacecraft 
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1 accelerometer results showing storm enhancements of oxygen density above auro- 

2 ral features by a factor of 3 to 8 [Bruinsma et al., 2006]. Despite the constraints cited 

3 above, it must be appreciated that this value is a parameter of the simple model 

4 used here. 

5 The available Poynting flux, evaluated at 4000 km, in conjunction with the iono- 

6 spheric conductance, determines the convection speed profile of the flux tube, which 

7 is then evaluated at the exobase. There, we assume the ring beam relaxes into a 

8 transverse power law, scaled by the local convection speed, as simulated and ob- 

9 served by Freja, in the form predicted by the Puhl et al., [1993] simulations. The 

10 fraction of the distribution that extends above the local escape speed, evaluated at 

1 1 the exobase, is taken to be the fraction of the local density that will escape. Thus, the 

12 partial density of the distribution that exceeds the escape speed (for the exponential 

13 velocity distribution, as an example] is: 


( 6 ] 


15 For a Lorentzian, the corresponding cumulative function is an arctangent, while for 

16 a Maxwellian, it is the complementary error function, erfc. These are used to com- 

17 pute the escaping fraction of the distribution, as summarized in Table 1. 
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The escape flux is estimated as the hemispheric flux, assuming that transverse en- 
ergy becomes parallel energy owing to magnetic flux tube divergence, and so is ob- 
tained as the product of the fractional density with the scale speed, taken to be v C on V : 
again for an exponential distribution, where Fii m = nVconv- 


F esc =n partial v cmv =F lim exp|- 


/ 1 

conv / 


( 7 ) 


The results of this integration of the escape flux are shown in Figure 4 for iono- 
spheric conductance of 8 S, in three cases, each with three sub-cases: 

First, we considered the result for escape of a generalized Lorentzian (power law) 
velocity distribution with scale speed equal to the convection speed at the exobase 
and with exponent matching the fit to the Freja data of Andre et al., [1994]. The 
downward shifted curve of similar shape indicates the result if the ambipolar poten- 
tial is ignored, while the nearly coincident dashed line indicates the same result if 
the centrifugal potential is ignored. 

Second, we considered an exponential velocity distribution with scale speed given 
by the convection speed at the exobase. Again, the downward shifted curve indicates 
the result if the ambipolar potential is ignored, while the nearby dashed curve indi- 
cates the result if the centrifugal potential is ignored. 
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Finally, we considered a ring beam that has relaxed fully to a Maxwellian with ther- 
mal speed given by the convection speed at the exobase, again with a parallel curves 
showing the results when ambipolar and then centrifugal potential is ignored. A plot 
of the Strangeway et al. [2005] data set and fit serves as background of Figure 4. 

The value of Poynting flux corresponding to D* = 0.3 is indicated on the plot to sug- 
gest a level at which the ring beam feature is appreciable, but as far as we know, no 
one has determined the threshold value of this parameter for wave growth. 

In Figure 5, we display the same results, in the case of the Lorentzian distribution, 
including centrifugal potential, for several different values of N p and corresponding 
ambipolar potentials. Insofar as the ambipolar potential is independent of the 
Poynting flux, the variations with N p explain some of the observed scatter in the es- 
cape flux data. A bold dashed curve cuts across the other traces as Poynting flux in- 
creases, obtained by requiring the precipitating electron density Np to track the 
Poynting flux according to the Strangeway et al., [2005] fits to the observations. This 
implements the reported fact that the two quantities are correlated, presumably re- 
sulting from common magnetospheric processes, with larger Poynting flux accom- 
panied by larger precipitating electron densities. Again, the value of Poynting flux 
corresponding to D* = 0.3 is indicated on the plot to show where the velocity distri- 
bution becomes perceptibly toroidal. 
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Figure 6 displays a plot of the dependence of escape flux on precipitating electron 
density, Np, for various assumed constant values of Poynting Flux. As in Figure 5, a 
bold dashed trace again represents the result of tying Np to Fp according to the 
Strangeway et al. [2005] fits. Clearly, independent behavior of these two parameters 
is capable of producing scatter in the results. Evidently, though the electron precipi- 
tation and ambipolar potential are both significant effects, they are not as powerful 
as the Poynting Flux in reproducing the observed large dynamic range of outflow 
flux. 

To summarize, little or no escape would be expected, at any but the extreme highest 
Poynting flux levels observed, after the distribution relaxes to a Maxwellian thermal 
speed corresponding to the local convective pick-up speed. Escape is appreciable at 
lower power levels if the velocity distribution at the exobase is an exponential with 
the local convection speed as scale speed. The escape response best mimics the ob- 
served (eyeball fit) empirical behavior for the Lorentzian power law distribution. 
When electron precipitation alone is assumed to drive the ambipolar potential, the 
result is a modest increase of escape flux at any Poynting flux level, an effect that 
nevertheless appears strong when coupled with a correlated Poynting Flux varia- 
tion. 


27 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


Moore & Khazanov - Mechanisms of Ionospheric Mass Escape 


Discussion 

Figure 7 is a schematic flow chart of the mechanisms described above, based loosely 
on the flow chart of Strangeway et al. [2005], but with the addition of an evaluation 
of the effects of ambipolar potential, centrifugal forcing and convective ion pick-up, 
and an allusion to neutral gas upwelling effects. The ambipolar electric field extends 
far above the exobase and is lumped into a total potential above the exobase. Cen- 
trifugal forcing also extends far above the exobase, lowering the height at which 
gravitational escape occurs. But while it may have strong effects on outflow speeds, 
we found that it only slightly reduces the escape speed for ionospheric ions at the 
exobase. The main effect of neutral gas upwelling is to raise the exobase, but we 
have not computed that effect here. 

Figure 7 also indicates the hypothetically important role of ion pick-up by convec- 
tion in producing unstable velocity distributions that drive waves to relax the ions 
into transversely accelerated power law distributions. The observed and modeled 
transverse velocity distribution at the exobase, acting in conjunction with the ambi- 
polar potential driven by precipitating electron thermal energy, implies increasing 
escape as ionospheric convection speed and electron precipitation increase, produc- 
ing wholesale escape of much of the velocity distribution for high values of the 
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1 Poynting flux that powers such convection, and-or of the precipitating electron den- 

2 sity that accompanies it. 

3 The results imply that the form of the ion velocity distribution maps directly into the 

4 observed dependence of escape on DC Poynting flux. A power law distribution gives 

5 the best agreement with the Strangeway et al. [2005] (power law] results, while a 

6 less extensive exponential distribution (in v] would suppress escape at lower levels 

7 of Poynting flux. Full thermalization of PUI ring distributions to a Maxwellian form 

8 would suppress outflow to well below the observed response at all power levels, 

9 even with substantial precipitating electron density, demonstrating the need for ac- 

1 0 tive wave particle ion acceleration. Variations of the actual velocity distribution 

1 1 produced are likely to be responsible for much of the scatter present in the observa- 

12 tions, though we have also shown that ambipolar potential variations are capable of 

13 contributing to the scatter. 

14 The curves in Fig. 5-6 give a reasonable eyeball fit to the FAST observations, sug- 

15 gesting that the power law distribution predicted by the simulations of Puhl et al. 

16 [1993] are realistic, as also observed in the ionosphere by the space shuttle Plasma 

17 Diagnostic Probe and Freja, at higher altitudes. Still it is clear that a more complete 

18 understanding of low energy ring beam relaxation via wave particle interactions is 

19 critically important to the proposed mechanism of plasma escape. The Puhl et al., 
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1 [1993] calculation is the only theoretical basis of which we are aware for ring beam 

2 relaxation into power law tails, as used here. As a semi-empirical parametric result, 

3 it is perhaps the weakest assumption in the proposed model, despite its successes. 

4 Further study of PUI ring beam relaxation is likely to yield significant results, possi- 

5 bly leading to a better understanding of the ubiquity of power law tails in space 

6 physics. If this paper motivates such study, it will have served its intended purpose. 

7 Four mechanisms are often invoked to explain ionospheric mass escape: i] heating 

8 of the F-region electron gas by addition of superthermal electrons, increasing the 

9 ambipolar potential, ii] Joule or frictional ion and neutral gas heating, raising scale 

10 heights and supply of plasma without escape, iii] resonant ion transverse accelera- 

1 1 tion (or ponderomotive forcing] by particle interactions with waves of unspecified 

12 source; and iv] centrifugal acceleration reducing the escape speed at the exobase. 

13 The results presented here suggest that convection of plasma relative to neutral gas 

14 is an important source of free energy (in toroidal PUI ring beams] for waves that ac- 

15 celerate ions near and below the exobase. Other wave sources are also likely to be 

16 important, but toroidal PUI distributions are a distinctive feature of ionospheric 

17 convection. The latter three mechanisms all respond to the speed of convection and 

18 thus to the amount of electromagnetic power (Poynting flux] from the linked mag- 

19 netospheric or solar wind plasmas. Mechanism iii] converts bulk motion into disor- 
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1 dered but non-thermal transverse ion energy that increases the number of ions that 

2 overcome gravity, as reduced by the ambipolar and centrifugal potentials. Mecha- 

3 nism iv) has been shown to operate at too high altitudes to be an significant con- 

4 tributor to escape fluxes, though it clearly accelerates outflows at higher altitudes 

5 and thus affects their ultimate destinations or escape from the magnetosphere. 

6 These "generalized Jeans escape" results suggest a theoretical basis for the empirical 

7 relationships of Strangeway [2005], identifying a new but important source for the 

8 waves that heat ionospheric ions, directly driven by convection and hence by DC 

9 Poynting flux. But this must be fully tested in future topside transport modeling. 

10 Also, a full understanding of the proposed process requires a solution of the quasi- 

11 linear relaxation of ring beams to determine its details and time scales and to fur- 

12 ther understand and validate the form of the velocity distribution that is formed. 

13 The hypothetical mechanism modeled here is motivated by the observed correlation 

14 of outflow with ionospheric convection as driven by DC Poynting flux. It is compati- 

15 ble by design with multiple features of auroral ion outflows, including their trans- 

16 verse heating in response to convection, response to electron precipitation, power 

17 law tails, and production of broadband ion-resonant waves. The results do not prove 

18 that there is no other source of the effective waves driven by Poynting flux, but they 
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1 do show that another source is unnecessary, given the operation of the mechanism 

2 as proposed. Two key tests of the proposed mechanism can be envisioned: 

3 If the asserted power law tails are rare at the exobase, or if their scale velocity corre- 

4 lates poorly with convection and Poynting flux, a different source of waves is im- 

5 plied. 

6 If another source of ion resonant wave energy is better correlated with DC Poynting 

7 flux or convection strength, it would refute this model. 

8 It has often been observed that powerful ionospheric escape events, referred to as 

9 ionospheric mass ejections [Moore et al., 1999], often contain appreciable amounts 

10 of the molecular ions N 2 + and NO + . Wilson and Craven [1998] showed that such 

11 events were associated with strong and prolonged (> 15 min] convection of the 

12 ionosphere. They argued that the heating required to lift molecular ions up to high 

13 altitudes might be provided by waves driven unstable by toroidal 0 + distributions. 

14 Such events are certainly increasing molecular scale heights, leading to formation of 

15 new molecular ions in sunlight or through charge exchange with the more abundant 

16 ion species. Thus it may also be profitable to consider the formation of toroidal pick- 

17 up distributions of the molecular ions. These should be more toroidal than those of 

18 0 + , since the molecular thermal speeds are smaller by a factor of ~2 than those of 0. 
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1 Thus it should be expected that enhanced molecular densities in the topside would 

2 lead directly to enhanced molecular ion upflow and outflow, assuming the mecha- 

3 nism suggested here is equally effective in generating power law molecular trans- 

4 verse velocity distributions. An assessment of molecular escape owing to ambipolar 

5 pick-up is beyond the scope of the present paper but would be a natural future out- 

6 growth of this work. 

7 Conclusions 

8 A model of ionospheric ion acceleration and outflow has been hypothesized on the 

9 basis of the following assumptions based on cited publications: 

10 i) Precipitating electrons are assumed to produce a superthermal secondary elec- 

11 tron population for which the ambipolar potential relation of Khazanov et al., [1997] 

12 applies. The fraction of the superthermal electron density above 50 eV is treated as 

13 a free parameter. 

14 ii) PUI ring beams form in the F region up to the exobase owing to convection of ions 

15 relative to neutrals, as described by St.-Maurice and Schunk [1979]. Ions are as- 

16 sumed weakly collisional but magnetized. The exobase height is an important pa- 

17 rameter. 
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1 iii) As D* becomes appreciable (>0.1-0. 3), PUI ring beams are assumed to generate 

2 broadband ion resonant waves that fill the depression at center of the ring beam, 

3 but also diffuse ions upward in energy (though always downward in phase space 

4 density), creating non-thermal power law velocity distributions with scale speed = 

5 V conv- 

6 Based on a simple "generalized Jeans escape" calculation, we find that escaping 0 + 

7 flux increases with Poynting flux (and associated convection) and with electron pre- 

8 cipitation in rough agreement with empirical scaling relationships of Strangeway et 

9 al., [2005]. The escape limit, Fij m , is treated as a parameter, the emphasis here being 

10 on the functional dependence on drivers rather than the absolute limiting flux. 
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i Table 1. Assumed velocity distributions 


Distribution type 

Definition (N = normalization) 

Cumulative 



( v esx < v < oo) 

Ring distribution 

/(Vj) = N ex p( - ( v / / 2 - (v ± - v co J 2 )lvl) 

v eS g = gravitational 

in plasma frame 


Vesc = g less centrifugal 

(v ! =Vconv)j 


Vesa = above less am- 

N = normalization 


bipolar 

Lorentzian (gen- 

m=N/(\ + (v L /v co j 2K ) 

2 .. 

l--arctan (v esx lv sd ) 

Jt 

eralized exponent, 



K) 



Exponential 

/(v ± )=iVexp(-v L /v conv ) 

exp(-V ra /V scl ) 

Thermalized 

/(v ± ) = Nex :p(-v 2 !v 2 conv ) 

erfc(Y e JVJ 
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i Figures 



4 electrons after Khazanov et al. [1997], digitized from their figure 3. 


5 
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Convection, Ion Escape Velocities 
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2 Figure 2. Altitude dependence of escape speed and convection speed for three val- 

3 ues of Poynting Flux at 4000 km altitude, at a typical value of ionospheric height in- 

4 tegrated conductivity (8 S). V esg is the normal gravitational escape speed. V eS c is the 

5 escape speed when centrifugal acceleration is taken into account. Circles mark 


6 points where the ring beam speed exceeds the local escape speed. 
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2 Figure 3. A representative example of the energy distribution of transversely accel- 

3 erated ions [After Andre et al. 1994], from the Freja mission at about 1700 km alti- 

4 tude over the auroral zone. Exponential and Lorentzian fits are shown as smooth 

5 curves. 
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Poynting Flux [DC] (mW/m 2 ) 

2 Figure 4. "Generalized Jeans escape" of 0 + escape flux vs. Poynting flux at 4000 km 

3 altitude, for three velocity distributions (Lorentzian, Exponential and Maxwellian 

4 ring beams), and for escape above three velocity thresholds: V eS g: gravitational only; 

5 Vesc: deducts centrifugal potential; V eS a deducts ambipolar potential. Data and fit 

6 from Strangeway et al., [2005] are included for comparison. Limiting flux is taken as 

7 2 x 10 10 cnr 2 s _1 . Exobase is taken at 500 km altitude, and ambipolar potential above 

8 500 km is 6.5 V corresponding to a precipitating electron density of 3/cm 3 4 5 6 7 8 9 and 

9 Np/Ntot = 0.3%. Poynting flux corresponding to D* = 0.3 is indicated for reference. 
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2 Figure 5. Similar to Figure 4, but with multiple traces reflecting Lorentzian velocity 

3 distributions, with the indicated values of precipitating electron density, Np. The 

4 trace for Np = 1 cnr 3 4 5 6 is the same as the Lorentzian trace on Figure 4. The bold 

5 dashed curve if for Np tracking PF per Strangeway et al. [2005], that is, for Np = 

6 (2.14 10 7 x PF 1265 /1.02 x 109)(i/2-2) 
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2 Figure 6. "Generalized Jeans escape" estimate of 0 + escape flux dependence upon 

3 precipitating electron density at 4000 km altitude, for multiple traces based on the 

4 indicated values of Poynting flux, assuming Lorentzian velocity distributions. Data 

5 and fit from Strangeway et al., [2005] are included for comparison. The limiting flux 

6 is taken as 2 x 10 10 cm 2 3 4 5 6 s -1 . The exobase is taken at 500 km altitude. 
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4 
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Figure 7. Schematic diagram of ionospheric outflow mechanisms illustrating the role 
of ion pick-up, ambipolar electric field, and centrifugal forcing in the outflow of 
ionospheric plasma. The flow chart is derived from that given by Strangeway et al. 
[2005]. 
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1 Figure Captions 

2 Figure 1. Ambipolar potential as a function of the fractional density of superthermal 

3 electrons after Khazanov et al. [1997], digitized from their figure 3. 

4 Figure 2. Altitude dependence of escape speed and convection speed for three val- 

5 ues of Poynting Flux at 4000 km altitude, at a typical value of ionospheric height in- 

6 tegrated conductivity [8 S). V eS cg is the normal gravitational escape speed. V eS cc is the 

7 escape speed when centrifugal acceleration is taken into account. Circles mark 

8 points where the ring beam speed exceeds the local escape speed. 

9 Figure 3. A representative example of the energy distribution of transversely accel- 

10 erated ions [After Andre et al. 1994], from the Freja mission, observed at about 1700 

11 km altitude over the auroral zone. Exponential and Lorentzian fits are shown. 

12 Figure 4. "Generalized Jeans escape" of 0 + escape flux vs. Poynting flux at 4000 km 

13 altitude, for three velocity distributions [Lorentzian, Exponential and Maxwellian 

14 ring beams), and for escape above three velocity thresholds: V eS g: gravitational only; 

15 Vesc: deducts centrifugal potential; V eS a deducts ambipolar potential. Data and fit 

16 from Strangeway et al., [2005] are included for comparison. Limiting flux is taken as 

17 2 x 10 10 cm 2 s _1 . Exobase is taken at 500 km altitude, and ambipolar potential above 
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1 500 km is 6.5 V corresponding to a precipitating electron density of 3/cm 3 and 

2 Np/Ntot = 0.3%. Poynting flux corresponding to D* = 0.3 is indicated for reference. 

3 

4 Figure 5. Similar to Figure 4, but with multiple traces reflecting Lorentzian velocity 

5 distributions, with the indicated values of precipitating electron density, Np. The 

6 trace for Np = 1 cm 3 is the same as the Lorentzian trace on Figure 4. The bold 

7 dashed curve if for Np tracking PF per Strangeway et al. [2005], that is, for Np = 

8 (2.14 10 7 x PFi-265/1.02 x lO 9 )! 1 ^) 

9 Figure 6. "Generalized Jeans escape" estimate of 0 + escape flux dependence upon 

10 precipitating electron density at 4000 km altitude, for multiple traces based on the 

1 1 indicated values of Poynting flux, assuming Lorentzian velocity distributions. Data 

12 and fit from Strangeway et al., [2005] are included for comparison. The limiting flux 

13 is taken as 2 x 10 10 cm 2 s _1 . The exobase is taken at 500 km altitude. 

14 Figure 7. Schematic diagram of ionospheric outflow mechanisms illustrating the role 

15 of ion pick-up, ambipolar electric field, and centrifugal forcing in the outflow of 

16 ionospheric plasma. The flow chart is derived from that given by Strangeway et al. 

17 [2005]. 
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